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G. Bellini,^ J. Benziger,^ D. Bick,'^ G. Bonfini,* D. Bravo, ^ M. Buizza Avanzini/ B. Caccianiga,^ L. Cadonati,^ 
F. Calaprice/ C. Carraro,* P. Cavalcante,'* A. Chavarria,^ D. D'Angelo/ S. Davini,® A. Derbin,^ A. Etenko/" 

K. Fomenko,"''' D. Franco/^ C. Galbiati/ S. Gazzana,^ C. Ghiano,'* M. Giammarchi/ M. Goeger-Neff/^ 

A. GorettiJ L. Grandi/ E. Guardincerri,® S. Hardy,^ Aldo lanni,'' Andrea lannij A. Kayunov,^ D. Korablev/^ 

G. Korga,"' Y. Koshio,"* D. Kryn/^ M. Laubenstein,^ L. Lewke/'^ E. Litvinovich/° B. LoerJ F. Lombardi,'' 

P. Lombardi/ L. Ludhova/ I. Machulin,i° S. Manecki,^ W. Maneschg," G. Manuzio,^ Q. Meindl/^ 

monti.i M. Misiaszek.i^''' D. Montanari.'*''^ P. Mosteiro7 V. Miiratova.^ L. Oberaiier.^^ 



r. Lombardi, L. ijuunuva, ±. iviauiiumi, o. iviaiici^Ri, vv . iviaiitBuiig, u. iviaiiuiiu, Vc^. iviciiiui, 

E. Meroni/ L. Miramonti,^ M. Misiaszek/^'^ D. Montanari,'*''' P. Mosteiro/ V. Muratova,^ L. Oberauer, 
M. Obolensky/2 F. Ortica/'^ K. Otis,*^ M. Pallavicini,^ L. Papp,^ L. Perasso/ S. Perasso,® A. Pocar,^ 

R.S. Raghavan,^ G. Ranucci/ A. Razeto,'^ A. Re/ P.A. Romani/'' A. SabelnikoVj^" R. Saldanha/ 

C. Salvo,^ S. Sch6nert/3 H. Simgen,!^ M. Skorokhvatov/° O. Smirnov," A. Sotnikov," S. Sukhotin/o 

Y. Suvorov,'' R. Tartaglia,"* G. Testera,^ D. Vignaud,^^ R.B. Vogelaar,^ F. von Feilitzsch,!^ j^ Winter/^ 

M. Wojcik/'^ A. Wright/ M. Wurm/ J. Xu/ O. Zaimidoroga," S. Zavatarelli/ and G. Zuzel^^ 

^ Dipartimento dt Fisica, Universita' degli Studi e INFN, 20133 Milano, Italy 

^ Chemical Engineering Department, Princeton University, Princeton, NJ 08544, USA 

■' Institut fur Experimentalphysik, Universitat, 22761 Hamburg, Germany 

^ INFN Laboraton Nazionali del Gran Sasso, SS 11 bis Km 18+910, 67010 Assergi (AQ), Italy 

^ Physics Department, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA 

''Physics Department, University of Massachusetts, Amherst, AM01003, USA 

'^ Physics Department, Princeton University, Princeton, NJ 08544, USA 

* Dipartimento di Fisica, Universita' e INFN, Genova 16 146, Italy 

^St. Petersburg Nuclear Physics Institute, 188350 Gatchina, Russia 

"' NRG Kurchatov Institute, 123182 Moscow, Russia 

^^ Joint Institute for Nuclear Research, 141980 Dubna, Russia 

^^ Laboratoire AstroParticule et Cosmologie, 75231 Pans cedex 13, France 

^'' Physik Department, Technische Universitaet Muenchen, 85747 Garching, Germany 

^"^ Max- Plank- Institut fuer Kernphysik, 69029 Heidelberg, Germany 

^^M. Smoluchowski Institute of Physics, Jagellonian University, 30059 Krakow, Poland 

^^Dipartimento di Chimica, Universita' e INFN, 06123 Perugia, Italy 

Borexino collaboration 

A search for 5.5-MeV solar axions produced in the p + d ^ ^He + A (5.5 MeV) reaction was 
performed using the Borexino detector. The Compton conversion of axions to photons, A+e — > e+7; 
the axio-electric effect, A.-\-e-\-Z ^- e-\-Z; the decay of axions into two photons, A — >■ 27; and inverse 
Primakoff conversion on nuclei, A + Z — >■ 7 + .^, are considered. Model independent limits on axion- 
electron [gAe], axion-photon [gA-t), and isovector axion-nucleon (gsAiv) couplings are obtained: 
\gAe X gsAivl < 5.5 X 10"" and I5A7 x gsAivl < 4.6 x 10""GeV"^ at rriA < 1 MeV (90% c.l.). 
These limits are 2-4 orders of magnitude stronger than those obtained in previous laboratory-based 
experiments using nuclear reactors and accelerators. 
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I. INTRODUCTION 

The axion hypothesis was introduced by Weinberg jl| 
and Wilczek [2], who showed that the solution to the 
problem of CP conservation in strong interactions, pro- 
posed earlier by Peccei and Quinn Q , should lead to the 
existence of a neutral pseudoscalar particle. The orig- 
inal WWPQ axion model produced specific predictions 
for the coupling constants between axions and photons 
[gA-^), electrons {gAe), and nucleons (^atv) which were 
soon disproved by experiments performed with reactors 
and accelerators, and by experiments with artificial ra- 
dioactive sources [J]. 

Two classes of new theoretical models, hadronic or 



KSVZ d, i] and GUT or DFSZ [7, "s^, describe "invis- 
ible" axions, which solve the CP problem in strong in- 
teractions and interact more weakly with matter. The 
scale of Peccei-Quinn symmetry violation {Ja) in both 
models is arbitrary and can be extended to the Planck 
mass mp « 10^^ GeV. The axion mass in these models 
is determined by the axion decay constant Ja'- 



niA ~ {Um^/fA){Vz/{l + z)), 



(1) 



where rn^ and /^ are, respectively, the mass and decay 
constant of the neutral tt meson and z = rriu/md is u and 
d quark-mass ratio. The equation ([T} can be rewritten as: 
TOA(eV) SS 6.0 X 10^//^ (GeV). Since the axion-hadron 
and axion-lepton interaction amplitudes are proportional 



to the axion mass, the interaction between axions and 
matter is suppressed. 

The effective couphng constants gA-y, gAe, and gAN 
are to a great extent model dependent. For example, the 
hadronic axion cannot interact directly with leptons, and 
the constant gAe exists only because of radiative correc- 
tions. Also, the constant gA-y can differ by more than 
two orders of magnitude from the values accepted in the 
KSVZ and DFSZ models Q. 

The results from present-day experiments are inter- 
preted within these two most popular axion models. The 
main experimental efforts are focused on searching for 
an axion with a mass in the range of 10~^ to 10~^ eV. 
This range is free of astrophysical and cosmological con- 
straints, and relic axions with such a mass are considered 
to be the most likely dark matter candidates. 

New solutions to the CP problem rely on the hypoth- 
esis of a world of mirror particles [id . Ill| and super- 
symmetry [IJ] . These models allow the existence of ax- 
ions with a mass of about 1 MeV, which are not precluded 
by laboratory experiments or astrophysical data. 

The purpose of this study is to search experimentally 
for solar axions with an energy of 5.5 MeV, produced in 
the p + d ^ ^He -|- A (5.49 MeV) reaction. The axion 
flux is thus proportional to the jPjP-n eutrino flux, which 
is known with a high accuracy [13l Il4l |. The range of 
axion masses under study has been extended to 5 MeV. 
The axion detection signatures exploited in this study are 
Compton axion to photon conversion, A-l-e — ?> e-l-7, and 
the axio-electric effect, A+e+Z — ^ e+Z. The amplitudes 
of these processes are defined by the gAe coupling. We 
also consider the potential signals from axion decay into 
two 7-quanta and from inverse Primakoff conversion on 
nuclei, A + Z — >■ 'y + Z. The amplitudes of these reactions 
depend on the axion- photon coupling gAy The signature 
of all these reactions is a 5.5 MeV peak. 

We have previously published a search for solar axions 
emitted in the 478 keV Ml-transition of ^Li using the 
Borexino counting test facility |15| . 

The results of laboratory searches for the axion as well 
as astrophysical and cosmological axion bounds can be 
found in [j]. 



II. 



THE FLUX OF 5.5 MEV AXIONS 



The Sun potentially represents an efficient and in- 
tense source of axions. One production mechanism is 
photon-axion conversion in the electromagnetic fields of 
the solar plasma. In addition, electrons could produce 
axions via Compton processes and bremsstrahlung. Fi- 
nally, monochromatic axions could be emitted in mag- 
netic transitions in nuclei, when low-lying levels are ther- 
mally excited by the high temperature of the Sun. 

Even the reactions of the pp-solar fusion chain and the 
CNO cycle can produce axions. The most intense fiux is 
expected from the formation of the "^He nucleus: 



According to the Standard Solar Model (SSM), 99.7% 
of all deuterium is produced from the fusion of two pro- 
tons, p + p— 7- d + e'^ + Ve, while the remaining 0.3% is 
due to the p + p + e^ -^ d-\- v^. reaction. The produced 
deuteron captures a proton with lifetime r = 6s. 

The expected solar axion flux can thus be expressed in 
terms of the pp-neutrino flux. The proportionality fac- 
tor between the axion and neutrino fluxes is determined 
by a dimensionless axion- nucleon coupling constant gAN , 
which consists of isoscalar g^AN and isovector g^AN com- 
ponents. The ratio between the probability of an Ml 
magnetic nuclear transition with axion production (oja) 
and photon production {ujy) can be expressed as |16|-p^: 



UJA_ 



1 



1 



27ra 1 + (52 



90AnPi + g^AN 



{fiQ - 0.5)/3i + /i3 - ?7i 



PA 

P-i 



(3) 

where p^ and pA are, respectively, the photon and ax- 
ion momenta; (5^ = E/M is the ratio between the 
probabilities of E and M transitions; a « 1/137 is 
the fine-structure constant; fJ-o — fJ^p + fJ^n ~ 0.88 and 
M3 = Mp ~ Mn ~ 4.71 are, respectively, the isoscalar and 
isovector nuclear magnetic moments; and /3i and rji are 
parameters dependent on the specific nuclear matrix el- 
ements. 

Within the hadronic axion model, the constants goAN 
and Q3AN can be written in terms of the axion mass 

Mm-- 

&JA 1 + Z + W 

= -4.03 X lQ-^{mA/leV), (4) 



ffsAAT = -xf CD + F) = 

2jA 1 + Z + W 

= -2.75 X lO-\mA/leY). 



(5) 



P 



^He + 7 (5.5 MeV). 



(2) 



where rriN ~ 939 MeV is the nucleon mass, and z = 
TTT-u/md = 0.56 and w = niu/rns = 0.029 are u, d and 
s quark-mass ratios. Axial-coupling parameters F and 
D are obtained from hyperon semi-leptonic decays with 
high precision: F=0.462 ± 0.011, D= 0.808 ± 0.006 [2(|. 
The parameter Sfs, characterizing the flavor singlet cou- 
pling is poorly constrained: (0.37 < Sjs < 0.53) and 
(0.15 < Sfs < 0.5) were found in [2l| and [22|, respec- 
tively. The values of the axion-nucleon couplings given 
in (|4]) and ([5]) are obtained assuming S'/s=0.5. The value 
of u- and d-quark-mass ratio z = 0.56 is generally ac- 
cepted for axion papers, but it could vary in the range 
(0.35 — 0.6) [4:]. These uncertainties in Sfs and z could 
cause the values of g^AN and g^AN to differ from Q and 
(O by factors of (0.4-1.3) and (0.9-1.9) times, respec- 
tively. 

The values of g^AN and g^AN in the DFSZ model de- 
pend on an additional unknown parameter, but have the 
same order of magnitude: they have (0.3 — 1.5) times the 
values of the corresponding constants for the hadronic 
axion. 



III. INTERACTION OF AXIONS WITH 
MATTER AND AXION DECAYS 
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FIG. 1. Ratio of the emission probabilities for axions and 7 
quanta {uja/<-Oj) in the p + d ^ "^He + 7 reaction (curve 1, 
left-hand scale); cross section of the Compton conversion and 
axio-electric effect for 5.5-MeV axions on carbon atoms for 
QAe = 1 (curve 2 and 3, right-hand scale). 



In the p+d — > ■^He-l-7 reaction, the Ml-type transition 
corresponds to the capture of a proton with zero orbital 
momentum. The probability, x, of proton capture from 
the S state at energies below 80 keV was measured in [2j] ; 
at a proton energy of- 1 keV, x = 0.55 (^^ = o.82). The 
proton capture from the S state corresponds is an isovec- 
tor transition, and the ratio uja/u^ from expression (3), 
therefore depends only on gsAN [13 • 
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The calculated values of the uja/^i ratio as a function 
of the axion mass are shown in Fig[T] The expected solar 
axion flux on the Earth's surface is then 

$A0 = $.pp(waM) = 3.23 X 10^" {g3AN)^{p A /Pi)H^) 

where ^i,vv = 6.0 x 10^°cm~^s^^ is the pp solar neutrino 
flux [l3J, uM ■ Using the relation between gsAN and itla 
given by ([5]) , the ^ao value appears to be proportional to 
2.44 X 10~^m\{pA/p^)^ , where tjia is given 
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A. Axion-electron interactions: Compton 
conversion and the axio-electric effect 

An axion can scatter an electron to produce a photon in 
the Compton- like process A -I- e — > 7 -I- e. The Compton 
differential cross section for electrons was calculated in 
|17| , |18| , |24| . The energy spectrum of the 7-quanta 
depends on the axion mass, while the spectra of electrons 
can be found from relation E^ = Ea — E-y. Here, Ea — 
5.49 MeV, which is the Q- value of the p{d, ■^He)7 reaction. 
The integral cross section corresponding to this mode is 
[H, [a, M: 



ccc = 



g\^a 2m'^{m + EA)y 
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Am^p^ 



y(m? +y) 



■In- 



Ea+Pa-, 



PAV 



Ea -pA 



(8) 



where pA and Ea are the momenta and the energy of 
the axion respectively and y — 2mEA + ti^. The di- 
mensionless coupling constant gAe is associated with the 
electron mass m, so that gAe — C'em/fA, where Ce is 
a model dependent factor of the order of unity. In the 
standard WWPQ axion model, the values /a=250 GeV 
and Ce=l are fixed and gAe ~ 2 x lO"*^. In the DFSZ 
axion models Ce = l/3cos^ /Sdfsz, where /3dfsz is an ar- 
bitrary angle. Assuming cos^ /3dfsz=l, the axion-electron 
coupling is gAe=2.8xlO~^^mA where itia is expressed in 
eV units. The hadronic axion has no tree-level couplings 
to the electron, but there is an induced axion-electron 
coupling at one-loop level [19i] : 



3na^m f E Ja 

gAe = „ f T7 In — 

27r/a V ^ "^ 



24 



31 



z + w A 
In — 



(9) 



where n is the number of generations, N and E are the 
model dependent coefficients of the color and electro- 
magnetic anomalies and A sal GeV is the cutoff at the 
QCD confinement scale. The interaction strength of the 
hadronic axion with the electron is suppressed by a factor 

The integral cross section ucc calculated for gAe — 1 
is shown in FiglT] For axions with fixed gAe (curve 2 in 
Fig. [T|), the phase space contribution to the cross section 
is approximately independent of tua for niA < 2 MeV 
and the integral cross section is: 



o'cc ~ gA 



2 X 4.3 X 10" 



25cm2. 



(10) 



in eV units. 



The other process associated with axion-electron cou- 
pling is the axio-electric effect A + e + Z —>■ e + Z (the 
analogue of the photo-electric effect). In this process the 
axion disappears and an electron is emitted from an atom 



with an energy equal to the energy of the absorbed ax- 
ion minus the electron binding energy Ei,. The cross 
section of the axio-electric effect on K-electrons where 
the axion energy Ea ^ Ei, was calculated in [24] and 
has a complex form; it is shown in Fig. [TJ The cross 
section has a Z^ dependence and for carbon atoms the 
cross section is aAe ~ g\^xl.3xl0~^^ cm^/electron for 
niA < 1 MeV. This value is more than 4 orders of magni- 
tude lower than for axion Compton conversion. However, 
thanks to the different energy dependence {ace ~ Ea, 
'^Ae "^ (-E'a)^^^^) and Z^ dependence, the axio-electric 
effect is a potential signature for axions with detectors 
having high Z active mass |25i [. 

For axions with a mass above 2m, the main decay mode 
is the decay into an electron-positron pair; A ^ e~^ + e~ . 
The lifetime of an axion in the intrinsic reference system 
has the form: 



Te+e 



8^/(g 



Am^). 



The probability of an axion to reach the Earth is 
P{mA,PA) = exp(-T//re+e-), 



(11) 



(12) 



where Tf is the time of flight in the reference system 
associated with the axion: 



^/ = 



LniA iTiA L 
cpA Ea Pc 



(13) 



Here L — l.b y. IQ^^ cm is the distance from the Earth 
to the Sun and /? — pa/ Ea is the axion velocity in terms 
of the speed of light. The condition Tf < O.lrg+e- (in 
this case, 90% of all axions reach the Earth) limits the 
sensitivity of solar axion experiments to gAe < (10^^^ _ 

10-") [21. 



Axion-photon interaction: axion decay and the 
inverse Primakoff conversion on nuclei 



If the axion mass is less than 2m, ^ — > e+ + e^ decay 
is forbidden, but the axion can decay into two 7 quanta. 
The probability of the decay, which depends on the axion- 
photon coupling constant and the axion mass, is given by 
the expression: 



T27 



647r 



9An 



(14) 



where 17^7 is an axion-photon coupling constant with di- 
mension of (energy)-^ which is presented as in p3|.|19l|: 



gAj 



2TTfA 
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2(4 + z + w) 
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0^77 (15) 



where E/N is a model dependent parameter of the or- 
der of unity. E/N = 8/3 in the DFSZ axion models 
(Ca77=0.74) and E/N = for the original KSVZ axion 
(Ca77=-1.92). 



The phase space for decay depends on m,\. For r27 
measured in seconds, gA-y in GeV"^, and niA in eV, one 
obtains: 



r27 - 1.3 X 10^ff^2^^3 ^ 3 5 ^ ^o24^-5^-2^_ (^g) 

The flux of axions reaching the detector is given by 

$A = exp(-T//T27)$Ao = exp{-Tfg\^m\/647r)<^Ao 

(17) 
where ^aq is the axion flux at the at the Earth in case 
there is no axion decay ^, T27 is defined by ([T4l [T6 |) . and 
Tf , given by (J13p is the time of flight in the axion frame 
of reference. Because of axion decay, the sensitivity of 
experiments using solar axions drops off for large values 
oigl^m\. 

The number of A 
Vis: 



27 decays in a detector of volume 



AT ^ VmA 

PcEAT2y 



(18) 



section is |1S|: 
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This leads, using the KSVZ model, to expected Borexino 
event rates like those shown in FigJ2]for different values 
of niA- As can be seen in the Figure, the expected event 
rate is peaked, with a drop-off at low ruA due to the 
lower axion decay rate in the detector, and a decrease 
at high m^A resulting from the reduced flux from axion 
decay in flight. The maximum A^7 corresponds to m,A 
= ((8/6)T27m^/(T//mA))^^®= 65 keV, where T27 and r/ 
are defined by ^ and (USD. 

Another process depending on gA-y coupling is the Pri- 
makoff photo-production on carbon nuclei A ~\- ^^C -^ 
7 + ^^C. The integral inverse Primakoff conversion cross 



(19) 



Because the cross section depends on the gAy coupling, 
the decrease in the axion fiux due to A — )■ 27 decays 
during their flight from the Sun should be taken into ac- 
count. The axion flux at the detector was calculated by 
the method described above. The atomic-screening cor- 
rections for ^^c were introduced following the method 
proposed in [1^ . The expected conversion rate in Borex- 
ino is shown in Fig 15] for different values of tua- 



C. Escape of axions from the Sun 

Axions could be captured within the Sun. The re- 
quirement that most axions escape the Sun thus limits 
the axion coupling strengths accessible to terrestrial ex- 
periments. Each of the 4 axion-matter interactions con- 
sidered in this paper contribute to these limits. 

The flux of 5.5 MeV axions on the Earth's surface is 
proportional to the pp-neutrino flux, as given in equa- 
tion ([7]), only when the axion lifetime exceeds the time 
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FIG. 2. The expected number of axion decays (1) and inverse 
Primakoff conversions on ^^C nuclei (2) in the 100 t per day for 
KSVZ axion model. Lines (3) and (4) show the corresponding 
curves under the assumption that axions do not decay during 
their flight from the Sun. 



of flight from the Sun and when the flux is not reduced 
as a result of axion absorption by solar matter. Axions 
produced at the center of the Sun cross a layer of ap- 
proximately 6.8 X 10"^^ electrons/cm^ in order to reach 
the Sun's surface. Axion loss due to Compton conver- 
sion into photons in the solar matter imposes an upper 
limit on qac after vifhich the sensitivity of terrestrial ex- 
periments using solar axions is reduced. The cross section 
of the Compton conversion reaction for 5.5-MeV axions 
depends weakly on the axion mass and can be written as 
f^CC ~ 9Ae ^ 4.3 X 10~^^cm^. For gAe values below 10"^, 
the axion flux is not substantially suppressed. 

The maximum cross section of the axio-electric effect 
on atoms is ctag ~ ffie^^l-^ ^ 10~^^cm^ (see Fig[T]for 
carbon). The abundance of heavy {Z > 50) elements 
in the Sun is ^ 10^^ in relation to hydrogen [2g|. If 
QAe < 10~^, the change in the axion flux does not exceed 
10%. 

The axion-photon interaction, as determined by the 
constant qaj, leads to the conversion of an axion into 
a photon in a field of nucleus. The cross section of the 



reaction is ape ~ 9ajZ^ x 1.8 x 10 ^^cm^. Taking into ac- 
count the density of ^H and '^He nuclei, the condition that 
axions efhciently escape the Sun imposes the constraint 



gA-y < 10 ^GeV ^. Constraint for the other elements 
are negligible due to their low concentration in the Sun. 

The axion-nucleon interaction leads to axion absorp- 
tion in a threshold reaction similar to photo-dissociation: 
A + Z ^ Zi + Z2- For axions with energy 5.5-MeV this 
can occur for only a few nuclei: ^^0,^'^C, and ^H. It was 
shown in [27!] that axio-dissociation cannot substantially 
reduce the axion flux for gAN < lO"'^. 

In all, the requirement that most axions escape the Sun 
sets these limits on the matter-axion couplings - gAe < 
10-6, gAj < lO-^GeV^^ and gAN < 10^^. 



IV. EXPERIMENTAL SET-UP AND 
MEASUREMENTS 

A. Brief description of Borexino 

Borexino is a real-time detector for solar neutrino spec- 
troscopy located at the Gran Sasso Underground Labo- 
ratory. Its main goal is to measure low energy solar neu- 
trinos via (^,e)-scattering in an ultra-pure liquid scintil- 
lator. At the same time, however, the extremely high 
radiopurity of the detector and its large mass allow it to 
be used to study other fundamental questions in particle 
physics and astrophysics. 

The main features of the Borexino detector and its 
components have been thoroughly described in |28|-|37|. 
Borexino is a scintillator detector with an active mass 
of 278 tons of pseudocumene (C9H12), doped with 1.5 
g/liter of PPO (C15H11NO). The scintillator is housed in 
a thin nylon vessel (inner vessel - IV) and is surrounded 
by two concentric pseudocumene buffers (323 and 567 
tons) doped with a small amount of light quencher 
(dimethyl phthalate - DMP) to reduce their scintilla- 
tion. The two buffers are separated by a second thin 
nylon membrane to prevent diffusion of radon coming 
from PMTs, light concentrators and SSS walls towards 
the scintillator. The scintillator and buffers are contained 
in a Stainless Steel Sphere (SSS) with diameter 13.7 m. 
The SSS is enclosed in an 18.0-m diameter, 16.9-m high 
domed Water Tank (WT), containing 2100 tons of ultra 
pure water as an additional shield against external 7's 
and neutrons. The scintillation light is detected by 2212 
8" PMTs uniformly distributed on the inner surface of 
the SSS. The WT is equipped with 208 additional PMTs 
that act as a Cerenkov muon detector (outer detector) to 
identify the residual muons crossing the detector. All the 
internal components of the detector were selected follow- 
ing stringent radiopurity criteria. 



B. Detector calibration. Energy and spatial 
resolutions. 



In Borexino, charged particles are detected by scintil- 
lation light induced by their interactions with the liquid 



scintillator. The energy of an event is related to the to- 
tal collected light by the PMTs. In a simple approach, 
the response of the detector is assumed to be linear with 
respect to the energy released in the scintillator. The 
coefficient linking the event energy and the total col- 
lected charge is called the light yield (or photo-electron 
yield). Deviations from linearity at low energies can be 
taken into account including the ionization deficit func- 
tion /(/cb, E), where fcs is the empirical Birks' constant. 
The detector energy and spatial resolution were stud- 
ied with radioactive sources placed at different positions 
inside the inner vessel. For relatively high energies ( >2 
MeV), which are of interest for 5.5 MeV axion studies, 
the energy calibration was performed with a ^^^Am-^Be 
neutron source. One can find a detailed description of 
the energy calibration in [32, [S^l ■ Deviations of the 7- 
peak positions from linearity was less than 30 keV over 
the whole energy range. The energy resolution scales 
approximately as {a/E) ~ (0.058 + 1.1 x 10-^E)/Ve 
where E is given in MeV units. The position of an event 
is determined using a photon time of flight reconstruction 
algorithm. The resolution of the event reconstruction, as 
measured using the ^^''Bi-^^^'Po /S — a decay sequence, is 
13±2 cm [33. 




C. Data selection 

The experimental energy spectrum from Borexino in 
the range (1.0-15) MeV, containing 737.8 live-days of 
data, is shown in Figl3l At energies below 3 MeV, the 
spectrum is dominated by 2.6 MeV 7's from the /3-decay 
of 208TI in the PMTs and in the SSS. 

The spectrum obtained by vetoing all muons and 
events within 2 ms after each muon is shown by curve 
2, Fig 131 Muons are rejected by the outer detector and 
by an additional cut on the mean time of the hits belong- 
ing to the cluster and on the time corresponding to the 
maximum density of hits. This cut rejects residual muons 
that were not tagged by the outer water Cherenkov de- 
tector and that interacted in the pseudocumene buffer 
regions (see [S^ for more details). 

To reduce the background due to short-lived isotopes 
(1.1 s ^B, 1.2 s ^Li, etc; see {33]) induced by muons, an 
additional 6.5 s veto is applied after each muon crossing 
the SSS (curve 3, FigE]). This cut induces 202.2 days of 
dead time that reduces the live-time to 535.6 days. 

In order to reject external background in the 5.5 MeV 
energy region a fiducial volume cut is applied. Curve 4 
of Fig|3] shows the effect of selecting a 100 ton fiducial 
volume (FV) by applying a cut R < 3.02 m. Addition- 
ally, a pulse shape-discrimination analysis based on the 
Gatti optimal filter [38| is performed: events with nega- 
tive Gatti variable corresponding to 7- and /3-like signals 
are selected (see [Slj for more details). This cut does not 
change the spectrum for energies higher than 4 MeV. 



FIG. 3. Energy spectra of the events and effect of the selec- 
tion cuts. From top to bottom: (1) raw spectrum; (2) with 
2 ms muon veto cut; (3) with events within 6.5 s of a muon 
crossing the SSS removed; (4) events inside FV 



D. Simulation of the Borexino response functions 

The Monte Carlo (MC) method has been used to sim- 
ulate the Borexino response S{E) to electrons and 7- 
quanta produced by axion interactions. The MC simu- 
lations are based on the GEANT4 code, taking into ac- 
count the effect of ionization quenching and non-linearity 
induced by the energy dependence on the event position. 
Uniformly distributed 7's were simulated inside the entire 
inner vessel, but only those which reconstructed within 
the FV were used in determining the response function. 
The MC candidate events were selected by the same cuts 
applied in the real data selection. 

The energy spectra of electrons and gammas from the 
axion Compton conversion were generated according to 
the differential cross section given in [13, [3, [2J] for 
different axion masses [l5|. The responses for the ax- 
ion decay into two 7 quanta were calculated taking into 
account the angular correlation between photons. The 
response functions for axion Compton conversion (elec- 
tron and 7-quanta with total energy of 5.5 MeV), for the 
axio-electric effect (electron with energy 5.5 MeV), ax- 
ion decay (two 7-quanta with energy 2.75 MeV in case of 
non-relativistic axions) and for Primakoff conversion (5.5 
MeV 7-quanta) are shown in FigHl The response func- 
tions are normalized to 1 axion interaction (decay) in the 



probabilities: 



L = 




W' 



-Nj 



{Nl^r- /A^fP! 



(21) 



where N!^^ and N![^^ are the expected (|20p and measured 
number of counts in the i-th bin of the spectrum, respec- 
tively. The dispersion of the peak (a) was fixed, while the 
position (£^o) was varied around Emc i 30 keV, to take 
into account the uncertainty in the energy scale. The 
others 4 parameters (a, b, c and S) were also free. The 
total number of the degrees of freedom in the range of 
3.2-8.4 MeV was 46. 

The fit results, corresponding to the maximum of L 
at 5=0 are shown in FiglSJ The value of modified 
X^ = El^r""-^*^ )V^*'' is X^= 44/46. Because of the 
low statistics, a Monte Carlo simulation of ([20)1 is used 
to find the probability of Xp ^ 44. The goodness-of-fit 
{p = 52%) shows that the background is well described by 
function (|20p . The upper limit on the number of counts 
in the peak was found using the imax('5') profile, where 
^max('S') is the maximal value of L for fixed S while all 
others parameters were free. The distribution of Lyaa,x{S) 
values obtained from the MC simulations for S > was 
used to determine confidence levels in iinax(»S'). The lim- 
its obtained on the number of events for different pro- 
cesses are shown in table 1. 



FIG. 4. Simulated responses to axion interactions in the 
Borexino IV: 1- axio-electric effect ( 5.49 MeV electrons), 
2- Compton a^cion to photon conversion (electrons and 7- 
quanta), 3- Primakoff conversion (5.49 MeV 7-quanta), 4- de- 
cay A — !■ 27. The inset shows the corresponding responses for 
events reconstructed within the FV. 



IV. The shift in the position of the total absorption peak 
for interactions involving 7's is caused by an ionization 
quenching effect. All response functions are fitted with 
Gaussians. 



V. RESULTS AND DISCUSSIONS 

A. Fitting procedure 

Figure [S] shows the observed Borexino energy spectrum 
in the (3.0 — 8.5) MeV range in which the axion peaks 
might appear. The spectrum is modeled with a sum of 
exponential and Gaussian functions, 



iV"X^)=a + fexe-"-^ + (S'/ 



2Tra) X e ^^ 



(20) 



where the position Emc (=5.49 MeV) and dispersion a 
(^0.15 McV) are taken from the MC response, S is the 
peak intensity and a, b and c are the parameters of the 
function describing the continuous background. 

The number of events in the axion peak S was calcu- 
lated using the maximum likelihood method. The likeli- 
hood function assumes the form of a product of Poisson 



TABLE I. The upper limits on the number of axions registered 
in Borexino FV (counts/536 days). CC - Compton axion to 
photon conversion, A + e ^- e + 'y; AE - axio-electric effect, 
A + e + Z — >■ e + Z; PC - Primakoff conversion on nuclei. 



A+^'C^-/- 



' C. The limits are given at 68(90)% c.l. 



reaction 


CC 


AE 


y4->27 


PC 


glim 


3.8 (6.9) 


3.4 (6.5) 


4.8 (8.4) 


3.8 (6.9) 



The limits obtained {S^jc - 0-013 c/(100 1 day) at 90% 
c.l.) are very low, e.g. ~ 10"* times lower than expected 
number of events from pp— neutrino (135 c/(100 t day)). 
The upper limits on the number of events with energy 
5.5 MeV constrain the product of axion flux ^a and the 
interaction cross section with electron, proton or carbon 
nucleus (TA-e,p.,c via 



•S'cvonts = ^AO'A-e,p,cNe.p.cTe < S 



lim 



(22) 



where N^^p^c is the number of electrons, protons and car- 
bon nuclei in the IV, T is the measurement time and e is 
the detection efficiency. The individual rate limits are: 



<^ACrA-e < 4.5 X lO^^^s^' 

'^ACTA-p < 2.5 X lO^^^s"^ 

^A<JA-c < 3.3 X lO^^^s"^ 



(23) 
(24) 
(25) 



These limits show very high sensitivity to a model- 
independent value $^cta- For comparison the standard 
solar neutrino capture rate is SNU = 10~^^s~"'^atom~^. 
A capture rate of solar neutrinos measured by Ga-Ge 
radiochemical detectors is about 70 SNU. 




E, MeV 




FIG. 5. The fitted Borexino spectrum in the (3.2 - 8.4) MeV 
range. Curve 3 is the detector response function for Compton 
axion-photon conversion at the 90% c.l. upper limit (5=6.9 
events) . 



B. Limits on qae and qan couplings 

The number of expected events due to Compton con- 
version in the FV of the detector are: 



Sec = '^ypp{^Al^-i)(^CcNeT£ 



(26) 



where ace is the Compton conversion cross sections, 
$^ = <^^pj,{ujAluJ.y) istheaxionfluxilTl), N^. = 9.17x10^^ 
is the number of electrons in the IV; T = 4.63 x 10^ s 
is the exposure time; and e — 0.358 is the detection effi- 
ciency obtained with MC simulations (FigH]). 

The axion flux ^a is proportional to the constant 
{93 an)'^, and the cross section ace is proportional to the 
constant g\^, according to expressions (O and ([8]). The 
Sec value depends, then, on the product of the axion- 
electron and axion-nucleon coupling constants: g^^ x 
{gsAN)^- According to Eqs. ([7]) and (ITUl) . and taking 
into account the approximate equality of the momenta of 
the axion and the 7-quantum {(pA/p-y)^ — 1 for ttia < 1 



FIG. 6. The limits on the qab coupling constant obtained by 
1- present work, 2 - present work for \gAe x gsAivl, 3- reactor 
[3^1 Efl] and solar experiments [13, [23|, 4- beam dump exper- 
iments [4ll . |43 | , 5- ortho-positronium decay EJ , 6- CoGeNT 
[4^ . 7- CDMS [ill, 8- solar axion luminosity ^TJ > 9-resonance 
absorption [4^, 10- read giant [40]. The excluded values are 
located above the corresponding lines. The relations between 
QAe and 771a for KSVZ- and DFSZ-models are shown also. 



excluded \gAe x gsAN] values is shown in FiglS] (line 2). 
At {pA/p-yY ~ 1 or niA < IMeV the hmit is: 



\gAe X 53Aw| < 5.5 X 10-^^ (90% c.l.) 



(28) 



The dependence of \gAe x gsAAfj on niA arises from the 
kinematic factor in equations ((6|) and ([8]); thus, these con- 
straints are completely model-independent and valid for 
any pseudoscalar particle. It's important to stress that 
the limits were obtained on the assumption that axions 
escape from the Sun and reach the Earth, which implies 
gAe < 10"*^ for TJiA < 2m and gAe < (10"" - lO^^^) jf 
ruA >2m ([H). 

Within the hadronic (KSVZ) axion model, gsAN and 
niA are related by expression ([5]), which can be used to 
obtain a constraint on the gAe constant, depending on 



MeV), the expected number of events can be written as: t^e axion mass (Fig[Sl line 1). For {pA/Pyf « 1 the 



Sec = 9Ae X glAN X 1.4 X lO^^^N.Te = 



limit on gAe and tua is: 



9Ae 



X glAN X 2.1 X 1025. (27) 

Using this relationship, the experimental 5^;™ can be 
used to constrain gAe x gsAN and mA- The range of 



\gAe X m^ I < 2.0 X 10^5 cV (90% c.l.), (29) 

where niA is given in eV units. For niA — 1 MeV, this 
constraint corresponds to gAe < 2.0 x 10"^^. Figure [Hj 



shows the constraints on gAe that were obtained in ex- 
periments with reactor, accelerator, and solar axions, as 
well as constraints from astrophysical arguments. 



C. Limits on qa-^ and qan couplings 



KSVZ model (TTSt . we exclude axions with mass in the 
(7.5 - 76) keV range (see Figl5]). Similar constraints 
can be obtained for DFSZ axions for specific values of 

COS^^dfez- 

The number of expected events due to inverse Pri- 
makoff conversion is: 



The analysis of ^ — >■ 27 decay and Primakoff photo- 
production is more complicated because axions can decay 
during their flight from the Sun. The exponential depen- 
dence of the axion flux on (7^7 and rriA, given by (IT71) . 
must be taken into account. 

The number of events detected in the FV due to axion 
decays into 2 7's within the IV are: 



Spc — ^Af^PcNcTepc 



(34) 



^2^ ~ N^Te2^ 



(30) 



where N^ is given by (jTS)) and £27 = 0.35 is the de- 
tection efficiency obtained by MC simulation. The re- 
leads to model-independent limits 



lation ^27 < S^X"2f 



on g^AN ^ dA-y ^^ axion mass. The expected value of ^27 
has a complex dependence on gA-y, gsAN and ttia given 
by equations p4)) - (|T8)) . 

In the assumption that /? « 1 the number of decays in 
the FV depends on glAN^ 9A-y and m^: 



N^ = 1.68 X 10 



-4„2 



5^7 X 93 AN X I^A, 



(31) 



where gA^y and niA are given in GeV and eV units, 
respectively. The limit derived from equation (j30p . at 
90% C.I., is 



IdA-y X gsANl 'xrn\< 3.3 x 10 " eV. 



(32) 



The dependence of 5*27 on gA-y and itia is obtained 
from ([5]), which gives the relationship between gsAN and 
niA in the KSVZ model. The relation ^27 < 5'^'!^27 "^" 
poses constraints on the range of gA-y and m^ values. 
The excluded region is inside contour la in FiglT] (90 % 
c.L). For higher values of g\~'m\ axions decay before 
they reach the detector, while for lower g\^'m\ ^^'^ prob- 
ability of axion decay inside the Borexino volume is too 
low. The limits on gA-f obtained by other experiments 
are also shown. 

The Borexino results exclude a large new region of 
axion-photon coupling constant (2 x lO"^**— 10~^)GeV^^ 
for the axion mass range (0.01 — 5) MeV. The Borex- 
ino limits are about 2-4 order of magnitude stronger 
than those obtained by laboratory-based experiments us- 
ing nuclear reactors and accelerators. Moreover, our ex- 
cluded region has begun to overlap the predicted regions 
from heavy axion models '1CHL3| . 

At niA < IMeV the constraint on gA-y and to a is given 
by 



\gA^\ X TO^ < 1.2 X 10"^ eV^ 



(33) 



So, e.g., niA—l MeV corresponds to gA-y < 1.2 x 
10~^^ GeV~^. Under the assumption that the axion- 
photon coupling gAy depends on axion mass as in the 



where apc is the Primakoff conversion cross sections; 
Nc is the number of carbon nuclei in the IV, and epc 
is the detection efficiency for 5.5 MeV 7's. The axion 
flux, $yi, is proportional to the constant g^AN^ ^^'^ ^^® 
cross section apc is proportional to the constant g^-y, ac- 
cording to equations ([7]) and (|19l) . As a result, the Spc 
value depends on the product of the axion-photon and 



axion-nucleon coupling constants: gj^ 



93AN- 



Under 



the assumption that $^ w <I>^o (true for 5^-y(GeV ^) x 
m^(eV) < 1.2 X 10^) one can obtain the limit: 

I5A7 X gsANl < 4.6 X 10-" GeV-i (90% c.L), (35) 

where again gA-y is in GeV^^ units. This limit is 25 times 
stronger than the one obtained by CAST [5^, which 
searches for conversion of 5.5 MeV axions in a laboratory 
magnetic field {\gAy x g^AN] < 1-1 x 10~^ at tua < leV). 
In the KSVZ model (O, the constraint on gAy and to a 
is given by the relation: 



\gAy\ xmA< 1.7 X 10" 



(36) 



For niA—l MeV, this corresponds to gAy < 1.7 x 
IQ-^ GeV^^. The region of excluded values of gAy and 
niA are shown in Fig[71 line lb; under the assumption 
that gA-y depends on ttia as in the KSVZ model p5)) we 
exclude axions with masses between (1.5 - 73) keV (see 
Figl2]). Our results from the inverse Primakoff process 
exclude a new region of gAy values at tua ~ 10 keV. 



VI. CONCLUSIONS 

A search for 5.5 MeV solar axions emitted in the 
p(d, '^IIe)A reaction has been performed with the Borex- 
ino detector. The Compton conversion of axions into 
photons, the decay of axions into two photons, and in- 
verse Primakoff conversion on nuclei were studied. The 
signature of all these reactions is a 5.5 MeV peak in the 
energy spectrum of Borexino. No statistically signifi- 
cant indications of axion interactions were found. New, 
model independent, upper limits on the axion coupling 
constants to electrons, photons and nucleons. 



-13 



and 



\gAe X gsANl < 5.5 X 10 



I5A7 X g3AN\ < 4.6 X 10-"GeV 



(37) 



(38) 



10 




10° 10° 10" 10"' 10° 10"= 



m^, MeV 



were obtained at ttia < IMeV and 90% c.l. 

Under the assumption that gsAN depends on m^ as 
in the KSVZ axion model, new 90% c.l. limits on axion- 
electron and axion-photon coupling as a function of axion 
mass were obtained: 



\gAe X tua] < 2.0 X 10"'' eV 



(39) 



and 



\gAj X m,A 



< 1.7 X 10 



'12 



(40) 



The new Borexino results exclude large regions of axion- 
electron and axion-photon coupling constants {gAe G 
(10-11 - 10-9) and gAj S (2 x W^^ - 10-^)GeV-i) 
for the axion mass range (0.01 — 5) MeV. 
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